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Abstract. In this paper, the asymptotic behaviour of the distribution tail of the station- 
ary waiting time W in the GI /GI /2 FCFS queue is studied. Under subexponentiahtype 
assumptions on the service time distribution, bounds and sharp asymptotics are given 
for the probabihty P{VF > x}. We also get asymptotics for the distribution tail of a 
stationary two-dimensional workload vector and of a stationary queue length. These 
asymptotics depend heavily on the trafhc load. 

Keywords: FCFS multi-server queue, stationary waiting time, large deviations, long 
tailed distribution, subexponential distribution. 

1. Introduction 

It is well known (see, for example, [151 [THl [1]) that in the stable single server 
first-come-first-served qneue GI /GI /I with typical interarrival time r and typical 
service time a the tail of stationary waiting time W is related to the service time 
distribution tail B{x) = P{(t > x} via the equivalence 

1 — 

P{VF>x}~- — / B{y)dy as x ^ 00, (1) 

provided the subexponentiality of the integrated tail distribution Bj defined by its 
tail 

B[{x) = min (1, 1^ B{y)dy), x>0. 

As usual we say that a distribution G on is subexponential (belongs to the 
class ^) if G * G{x) ~ 2G{x) as x — )• 00. The converse assertion is also true, 
that is, the equivalence ([T|) implies the subexponentiality of Bj, see [151 Theorem 
1] for the case of Poisson arrival stream and [141 Theorem 1] for the general case. 

In this paper we consider the GI /GI / s FCFS queue which goes back to Kiefer 
and Wolfowitz [13]. We have s identical servers, i.i.d. interarrival times {t„} with 
finite mean a = Eri, and i.i.d. service times {cr„} with finite mean b = Eo"i. The 
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sequences {r„} and {(T„} are mutually independent. The system is assumed to be 
stable, i.e., p = b/a £ (0,s). We are interested in the asymptotic tail behaviour 
of the stationary waiting time distribution P{VF > x} as x — t- oo. 

It was realized recently (see, for example, existence results for moments in 
[T6] . [T7] : an asymptotic hypothesis in [19]; asymptotic results for fluid queues 
fed by heavy-tailed on-off flows in [5]) that the heaviness of the stationary waiting 
time tail depends substantially on the load p in the system. More precisely, it 
depends on p via the value of /c G {0, 1, . . . , s — 1} for which k < p < k + 1. In 
particular, Whitt conjectured that 



"where 7 and ry are positive constants (as functions of x)" [sic, [19]] . In the present 
paper we show that, in general, the tail behaviour of W is more complicated. 

Let R{w) = {Ri{w), . . . , Rs{w)) be the operator on which orders the 
coordinates of u; G in ascending order, i.e., Ri{w) < • • • < Rs{w). Then the 
residual work vector Wn = (Wni, • • • , Wns) which the nth customer observes just 
upon its arrival satisfies the celebrated Kiefer-Wolfowitz recursion: Wi = i ■ 0, 

Wn+l = R{{Wnl +<Tn- Tn+l) + , {Wn2 - Tn+l) + , • • • , (Wns " ^^+1) + ) 
= R{Wn + eiGn - iT„+i) + , 

here ei = (1,0, ...,0), i = (1,...,1) and = (max(0, . . . , max(0, w^)). 
The value of Wni is the delay which customer n experiences. In particular, the 
stationary waiting time is a weak limit for Wni- 

The process Wn is a Markov chain in R*. It is well known that, for general 
multi-dimensional Markov chains, large deviation problems are very difficult to 
solve even for stationary distributions. Usually they can be solved in low dimen- 
sions only, 2 or 3 at most, see [12^ [4], Almost all known results are derived for 
so-called Cramer case which corresponds to light-tailed distributions of jumps. 
In the heavy-tailed case almost nothing is known for general multi-dimensional 
Markov chains. 

The process Wn presents a particular but very important example of a Markov 
chain in R*, even if we are interested in the first component Wni- As follows from 
our analysis, the case s = 2 can be treated in detail. The stability condition for 
this particular case is 6 < 2a. One of the following cases can occur: 

(i) the maximal stability case when b < a; 

(ii) the intermediate case when b = a; 

(iii) the minimal stability case when b G (a, 2a). 

We find the exact asymptotics for P{T1/^ > x} in the maximal and minimal 
cases. We also describe the most probable way for the occurrence of large de- 
viations. In the intermediate case, we only provide upper and lower bounds. 
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Then we study the asymptotics for the tail of the distribution of a stationary 
two-dimensional workload vector and give comments on the tail asymptotics of 
the stationary queue length. 

For s > 2, the stability condition is 6 < sa. We hope that, for s > 2, direct 
modifications of our arguments may lead to exact asymptotics in two particular 
cases when either b < a (the maximal stability) or 6 G {{s — l)a, sa) (the minimal 
stability). However, one has to overcome many extra technicalities for that. 
Insofar as the case b £ [a, (s — l)a] is concerned, we are extremely sceptical on 
the possibility to get any sharp tail asymptotics in explicit form. 

For the two-server queue, in the maximal stability case, we prove the following: 



Theorem 1. Let s = 2 and b < a. When the integrated tail distribution Bj 
is subexponential, the tail of the stationary waiting time satisfies the asymptotic 
relation, as x ^ oo, 

P{W > x} ~ — {Bi{x)f + b / Bj{x + ya)B{x + y{a - b))dy 

a[2a — o) L Jq 

The proof follows by combining the lower bound given in Theorem [3] (Section 
[3]) and the upper bound given in Theorem H] (Section H]). Simpler lower and upper 
bounds for P{iy > x} are given in the following 



Corollary 1. Under the conditions of Theorem^ 

2a + b , P{W>x} ^ P{W>x} 1 

< lim mi ,— , , , ^ < lim sup — . . „ < 



2a2(2a-6) " ^^oo {Bi{x))^ 



>oo {Bi{x)Y ~ 2a{a - b) 



In our opinion, in Theorem [T] it is possible to obtain a compact expression 
for the tail asymptotics of Y*{W > x) only in the regularly varying case. A 
distribution G (or its tail G) is regularly varying at infinity with index 7 > 
(belongs to the class ^^Y), if G{x) > for all x and, for any fixed c > 0, 
G{cx)/G{x) — >• c~'^ as X — )■ 00. 

Corollary 2. Let b < a and the tail distribution B of service time be regularly 
varying with index 7 > 1. Then 

P{W > x} ~ c'(Bi{x) f, 



where 



1 + 



dz 



a{2a-b)l 7 - 1 7o {1 + zap-'^il + z{a - b)p 
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Recall definitions of a number of classes of heavy-tailed distributions. A 
distribution G is long-tailed (belongs to the class ^) if G{x) > for all x and, 
for any fixed t, 

G{x + 1) 

— =— > \ as X — )• oo. 

G{x) 

A distribution G belongs to the class J' of intermediate regularly varying 
distributions if G{x) > for all x and 

limliminf ^ = 1. 

dl x^oo G{x) 

Clearly, C .J^^.Y. 

In the minimal stability case, we prove the following 



Theorem 2. Let s = 2 and a < b < 2a, B £ y and Bj G J^^Y. Then 

P{W > x} ~ — —-rBil-r— — x] as x ^ oo. 

2a — b ^b — a / 

The proof is given in Section [7] and is based on the lower and upper bounds 
stated in Theorems [5] and [6] respectively. 

One can provide simple sufficient conditions for B € ^ and Bj S J" S^'V . Let 
Si be the class of all distributions G on such that G{x) > for all x and 
lim inf ^^oo G(2x) /G(x) > 0. Then the fohowing are known: (i) M'f C J ^Y C 
^) C y\ (ii) if B G ^ has a finite first moment, then Bj G J (see e.g. 
[6]). Therefore, if G ^ ™d has a finite first moment, then B satisfies the 
conditions of Theorem [2j Note that the converse is not true, in general: there 
exists a distribution B G ^ with a finite first moment such that Bj G ^ M"V , 
hut B ^yf]S (see Example 2 in HI Section 6]). 

The paper is organized as follows. Section [2] contains some auxiliary re- 
sults. In Section [3l we formulate and prove a result concerning a lower bound 
for P{VF > x} in the maximal stability case. The corresponding upper bound 
is given in Section HI Sections [5l El and [7] deal, respectively, with lower bounds, 
upper bounds, and asymptotics for P{M^ > x} in the minimal stability case. In 
Section [HI we prove further results related to the joint distribution of a stationary 
workload vector. Comments on the asymptotics for a stationary queue length 
distribution may be found in Section [H 

A number of upper and lower bounds for P{VF > x} in s-server queue are 
proposed in Remarks O [Sj HJ and [H 
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2. Preliminaries 

[21.1. Reduction to deterministic input stream case in assertions asso- 
ciated with upper bounds. Consider a general GljGlj s queue. Take any 
a' G {h/s^a). Consider an auxiliary D/GI/s system with the same service times 
{(7,1 } and deterministic inter arrival times = a': W[ = and 

<+i = i?« + eia„-m')+. 

Let W be a stationary waiting time in this auxiliary system. 

Lemma 1. If'P{W' > x} < G{x) for some long-tailed distribution G, then 

hmsup — =— < 1. 

x-s-oo G[X) 

Proof. Denote = o,' — Tn- Put Mq = and, for n >1, 

Mn = max{0, in + (,n-l, Cn-\ h6} 

= max(0, e„ + M„_i) = (e„ + M„_i)+. 

First, we use induction to prove the inequality 

W„ < W;; + iMn a.s. (2) 

Indeed, for n = 1 we have < + iMi . Assume the inequality is proved for some 
n; we prove it for n + 1. Indeed, 

Wn+i = RiWn + eio-„ - ir„+i)+ 

< R{W;, + iMn + eian - ir„+i)+ 

= RiW^ + eian - ia' + i{Mn + in+i))^ ■ 

Since (u + w)+ < n+ + 

< + eic7„ - + i{Mn + = W;;+i + iMn+i, 

and the proof of ([2]) is complete. 

Let M be the weak limit for M„ which exists due to E,^i = a' — a < and 
Strong Law of Large Numbers. The following stochastic equality holds: 

M =st max{0, 6, 6 + 6, • • • , 6 + • • • + ^n, • • •}• 

Since the random variable is bounded from above (by a'), there exists /3 > 
such that Ee'^^^ = 1. Then by Cramer estimate (see, for example, [U Section 5]), 
for any x, 

P{M > x} < e"^^. (3) 
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The inequality ([2]) implies that W <st W'+M, where W and M are independent. 
Let a random variable r] have distribution G and be independent of M. Since 
V >st W , we have W <st r] + M. Therefore, for any h > 0, 

rx — h 

P{W >x} < P{M >x- y}P{r] G dy} + P{r/ > x - /i} 

JO 

< r e-''(^-?^)G((iy) + G(x-/i), 



by ([3]). Integrating by parts yields 

Jo 

< e"'^^' + /3 r G(x - y)e-^ydy. 
Jh 

The distribution G is long-tailed, thus, for any e > there exists such that 

G(x - 1) < G(x)e" 
for any x > x(e). Hence, there exists c(e) < oo such that 

G{x -y) < c{e)G{x)e'y 
for any x > x(e). Take e = /3/2. Then 



G(x - y)e-^ydy < c{e)G{x) / e-^^^/^^y < =^G{x)e 



/3/2 
Hence, 

P{W > x} < e"^^ + 2c(e)G(x)e-^'^/2 ^ _ 
Taking into account also that G{x — /i) ~ G{x) for any fixed h > 0, we obtain 
limsup^i^-4^ < 2c(e)e-/5'^/2 + 1. 

Letting /i — )• oo yields the conclusion of the Lemma. 

[2].2. Reduction to deterministic input stream case in assertions as- 
sociated with lower bounds. Take any a' > a. As in the previous subsection, 
consider an auxiliary D/GI/s system with the same service times {(Tn} and de- 
terministic interarrival times t!^ = a' . Let W' be a stationary waiting time in this 
auxiliary system. 
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Lemma 2. IfF{W' > x} > G{x) for some long-tailed distribution G, then 

liminfP'f >/' > 1. 

^^oo G(X) 

Proof. Put = rn- a', Mq = and 

M„ = max{0, ^n, + Cn-l, • • • , + • • • + 6} = (M„_i + ^n)^ . 

For any n > 1, the following inequality holds: 

Vr„ > < - iM„. (4) 

Indeed, by induction arguments, 

Wn+i = R{Wn + eian - ir„+i)+ 

> R{W;, - iMn + eian - ir„+i)+ 

= R{Wl^ + eian - ia' - i{Mn + ^n+i))^ ■ 

Since (u — v)~^ > — , 

Wn+i > R{W;^ + eian - ia')+ - i(M„ + = W;,^^ - iM„+i, 

and the proof of (jlj) is complete. 

Let M be the weak limit for M„ which exists due to E^i = a — a' < and the 
Strong Law of Large Numbers. The inequality ^ implies that W >st W — M 
where W and M are independent. Therefore, for any h > 0, 

P{W >x} > P{W' > x + h}P{M <h}> G{x + h)P{M < h}. 

The distribution G is long-tailed, thus G{x + /i) ~ G{x) for any fixed h > and 

liminf ^ > P{M < h}. 

Letting /i — )• oo, we obtain the desired estimate from below. 

[2l3. Adapted versions of the Law of Large Numbers. It is well known 
that obtaining lower bounds for systems under assumptions of heavy tails usually 
requires some variant of the Law of Large Numbers. Here we provide such a tool 
for the two-server queue. 

Lemma 3. Let {in^fln), n = 1, 2, . . . , be independent identically distributed 
pairs of random variables. Let the two-dimensional Markov chain Vn = {Vni, Ki2)) 
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n = 1, 2, . . . , be defined in the following way: Vi has an arbitrary distribution 
and 



Vr 



n+1 



Vn + {in,rin), if Vnl < Vn2, 
Vn + (f?n,Cn), if Vnl > K2- 



If^rji < E^i, then the following convergence in probability holds: 



n 



as n ^ oo. 



Proof. Since K+i,i + K+i,2 = Ki + + U + Vn, by the Law of Large 
Numbers 

> E^i + Er/i as n — )• 00. (5) 

n 

Define a Markov chain C/„ = Vn2 — Vni- If Un > 0, then — Un = r]n — ^n, 
while if [/„ < 0, then Un+i — Un = £,n — fin = —{Vn — S,n), so, Un is the oscihating 
random walk. Since E^i > Er/i, the mean drift of the chain Un is negative on 
the positive half-line and is positive on the negative half-line. Therefore, for any 
sufficiently large A, the set [—A, A] is positive recurrent for this Markov chain. 
In particular, the distributions of Un are tight. Hence, Un/n — )• in probability 
as n — )• 00. Together with ([5]), it implies the desired assertion of Lemma. 

The classical Law of Large Numbers and Lemma [3] imply the following 

Corollary 3. Let Brji < E^i < and e > 0. Then 

P{Vni > 0, 142 > 0| 1/1 = (^1,^2)} ^ 1 
as N ^ 00 uniformly in n > N and in {vi,V2) on the set 

|ui,W2 > n(|E^i| +e), vi + V2> n(|E^i| + |E?7i| +e)}. 
Corollary 4. Let Ern < E^i < and e > 0. Then 

P{Vnl > 0, K2 > 0| Fi = {vi,V2)} ^ 

as N ^ 00 uniformly in n > N and in {vi,V2) on the complementary set 
{vi > n(|E^i| - e), V2 > n(|E^i| - e), vi + V2 > n(|E^i| + |Er/i| - e)}. 
Corollary 5. Let Brji < 0, E^i > 0, Er/i + E^i < and e > 0. Then 

P{Kl > X,Vn2> x\Vi = {vi,V2)} 1 

as X, N ^ 00 uniformly in n > N and in (^1,^2) on the set 
> x-n(E^i-£), V2 > 2a; + n(|E^i + E7?i| 
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3. The maximal stability case: a lower bound 

Theorem 3. Assume b S (0, a). Let the integrated service time distribution 
Bj be long-tailed. Then the tail of the stationary waiting time W admits the 
following estimate from below: as x ^ oo, 



P{W>x}> ' + ^(') 



a{2a - b) 



{Bi{x)f + b / Bi{x + ya)B{x + y{a - b))dy 
Jo 



(6) 



Remark 1. From ([6]), one can get the lower bound in Corollary [TJ Namely, 
replace B{x + y{a — b)) by a smaller term B(x + ya) in the integral in the RHS of 
([6]). Then the new integral is equal to 6(i3/(x))^/2a, and the lower bound follows 
since 

1 / 6 \ _ 2a + 6 
a(2a - 6) V ^ 2^) ~ 2a'^{2a-b)' 

R e m a r k 2. By use of Strong Law of Large Numbers, one can get the following 
result for s-server queue, s > 2. If 6 < a, then there exists a constant K = 
K{a, b, s) such that 

P{W >x}>{K + o{l))(Bi{x)y. 

We start with some auxiliary results. The proof of the theorem is given in 
subsection [314. 

[3].l. An integral equality. 

Lemma 4. Let f[y) be an integrable function. Put fi{y) = f{z)dz. 
Then, for any positive a and j3, a > j3, 

coo roo 



roc roc 

J = 1 / f{ay+Pz)f{(3y+az)dydz 
Jo Jo 

ifm? 2/3 



Proof. Put u = ay+/3z and v = f3y+az. Then 



fi{au)f{(3u)du. 



]^ ^oo rau/fi 

^ = 2 M / /(^)^^ / f{v)dv 

a^-p^ Jo Jpu/a 

f{u)fiiPu/a)du I fiu)fi{au/P)du 



= 2 fl2 / f{au)fi{l3u)du / f{f3u)fi{au)du. 

a'^-p'^ Jo a'^-p'^ Jo 

Integration by parts yields 

fi{Mf{au)du = iM^ll! _ ^ r fj{au)f{(3u)du. 
a a Jo 
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By substituting this equality into the previous one, we arrive at the conclusion 
of the Lemma. 

[3].2. Some calculations vi^ith two big service times. Fix e > and 

put b' = h — £. For k and I, k < I < n, define the events Anki and Cnki by the 
equalities 

A-nki = {o'fc > X + [l - k)a + {n - l){a - b'), ai > x + (n - l){a - b'), 

ak + cri > 2x + {I - k)a + (n - l){2a - 

and 



n 

Cnki= f] {aj <x + {n-j){a-b')}. 



Note that the events Anki H C„fc; are disjoint for different pairs {k,l). Due to the 
existence of Eo", uniformly in n > 1 and k < I < n, 

oo 

as X oo. (7) 

j=Q 

Lemma 5. Assume b E (0, a). Let the integrated tail distribution Bj be long- 
tailed. Then, for any fixed N > 1 and for any e > 0, as x ^ oo, 

n-N 1 r i-oo 

{Bi{x)f + b' / Bj{x + ya)B{x + y{a - b'))dy 
Jo 



lim V P{Anki} ^ 



n->oo 

fc,i=i 
k<l 



a{2a - b') 



Proof. Put 

A'ki = {(Ji > X + ka + l{a—b'), a2> x + l{a—b'), ai + a2 > 2x + ka + /(2a— 6')}, 
so that P{Anki] = P{^;_fe,„_J and 

n—N n—1 n—l—1 oo oo 

„ii- E nAnki} = „iim E E nA'ki} = E E n^'ki}. (s) 

fc,i=i i=7V fc=l l=Nk=l 

k<l 

Consider also the events 

A{y,z) = {cji > x+ya+z{a—b'), a2 > x+z{a—b'), ai+a2 > 2x+ya+z{2a—b')}, 

which satisfy A(k, I) = A'^^. Since the probability Y'{A{y, z)} is non-increasing in 
y and z, we have the inequalities 

roo roo 



/_ = / / P{A{y, z)}dydz < P{^h} 

roo roo 

< / / P{A{y,z)}dydz = I+. (9) 
Jo Jo 
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The values of integrals /_ and I-^- are close to each other in the following sense: 



r'TV poo poo pi 



rl\ roo roo ri 

< / P{A{y,z)}dydz+ / P{A{y, z)}dydz 
Jo Jo Jo Jo 

< A^P{cJ2 > x} P{ai > x + ya}dy + P{cri > x} / P{cji > x + z{a-b')}dz. 

Jo Jo 

Recall that the distribution Bj{x) is long tailed, which is equivalent to B{x) = 
o{Bi{x)). Therefore, as x — )• oo, 

/+-/_< ^i±^B{x)Bi{x) = o{(Bj{x)f). 

Now it follows from ([9|) that, as x — )• oo, 

.z)}dydz + o{(Bi{x))^). (10) 



l=Nk=l 



Further, 



p{A{y,z)} 

= B{x + ya + za)B{x + z{a — h')) 

+ p|x + ya + z{a — b') < ai < x + ya + za, a2 > x + z{a — b'), 

+ o"2 > 2x + ya + z{2a — 6')| 

= 'B{x + ya + za)B{x + z{a - b')) 

+ p|x + ya + z{a—b') < ai < x + ya + za, ai + a2 > 2x + ya + z(2a— 6')| 
= B{x + ya + za)B{x + z{a - b')) + Q{y, z), 

since the event {fii < x + ya + za, o"! + o"2 > 2x + ya + z{2a — 6')} implies the 
event {a2 > x + z{a' — b)}. Consequently integrating over y and z, we obtain 

roo roo roc 

/ / B{x+ya+za)B{x+z{a—b'))dydz = — Bj{x+za)B{x+z{a—b'))dz. 
Jo Jo a Jo 

By the total probability formula, 

Q{y,z) = / P{ai £ X + ya + z{a-b') + dt}P{a2 > X + za - t} 
Jo 

rzb' _ 

= / B{x + za -t)B{x + ya + z{a-b') + dt). 
Jo 
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The integration against y leads to the equahties 



/ Q{y,z)dy = - B{x + za - t)Bi{x + z{a - b') + dt) 
Jo a Jo 

1 r^'>' — — 

= - / B(x + za- t)B{x + z(a - b') + t)dt 
a Jo 

b' — — 

= - / B(x + za- tb')B(x + z(a - b') + tb')dt. 
a Jo 



Integrating against z, we obtain: 



/ / Q{y,z)dydz = - / B{x + za - tb')B{x + z{a - b') + tb')dtdz 

Jo Jo a Jo Jo 

= - / / B{x + za- tb')B{x + z{a - b') + tb')dzdt 

a Jo Jt 

y roo roo 

= -/ / B{x + za + t{a-b'))B{x + z{a-b') + ta)dzdt. 

a Jo Jo 



By Lemma m with f{y) = B{x + y), a = a, and (3 = a — b', the latter integral is 
equal to 

Putting everything together into (fTOj) . we obtain the following equivalence, as 
X — )• oo: 

oo oo 



I)' roo 

"^ a(2a - b') Jo + + ~ 

which due to ^ completes the proof of Lemma. 

[3l3. Proof of Theorem [Si If Bi{x) is long-tailed, then the function in x 

^oo 

{Bi{x)f + b / Bi{x + ya)B{x + y{a - b))dy 
Jo 

is long-tailed as well. Indeed, for any fixed t, we have, as x — )• oo, 

roo roo 

/ Bi{x+t+ya)B{x+t+y{a-b))dy ~ / B i{x+ya)B{x+t+y{a-b))dy. 
Jo Jo 

Integrating by parts we get the equality for RHS integral 



a — b 

1 



oo 



^ Bi{x+ya)Bi{x+t+y{a-b)) - / B{x+ya)Bi{x+t+y{a-b))dy 

JO 



a — b 



^oo 

{Bj{x)f - / B{x+ya)Bi{x+y{a-b))dy 
Jo 

Bi{x+ya)B{x+y{a-b))dy. 
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So, we can apply Lemma [21 and it is sufficient to prove the lower bound of 
Theorem [3] for the queueing system D/GI/2 with deterministic input stream. 
Let the interarrival times r„ be deterministic, i.e., Tn = a. Then the event A^ki 
implies the event 

> x + {l-k)a + {n-l){a-b') - a, VFj+i^i > x + {n-l){a-h') - a, 

Wi+i^i + Wk+1,2 >2x + {l-k)a + (n-/)(2a-6') - 2a}, 

which implies 

W^m,i >x + {n-l){a-b')-a, Wi+i,i+Wi+i,2 > 2x + {n-l){2a-b')-2a} . 
Thus, by Corollary [3] (with ^ = o" — r and r] = — r), there exists N such that 

P{Wn > X I Ankl} > 1 - e (11) 

for any n > N and k < I < n — N . 

Taking into account that the events Anki H Cnki are disjoint for distinct pairs 
{k,l), we obtain the following estimates: 

n-N n-N 

^{Wn >X]>Y. > ^.AnkuCnkl] 
k=l l=k+l 

n-N n-N n-N n-N 

^ E E 'P{Wn>X,A^kl]-Y. E ^{^nkuCnkl}. 
k=l l=k+l k=l l=k+l 

Since the events Anki and Cnki are independent, 

n-N n-N n-N n-N 

^{Wn>x}> J2 E P{Wn>X,Ar,kl}-SUpP{Cnkl}Yl E ^{Ankl} 
k=l l=k+l k=l l=k+l 

n-N n-N n-N n-N 

= > ^ I Ankl}P{Ankl} - o(l) E ^{^nkl} 

k=l l=k+l k=l l=k+l 

as rr — )• oo uniformly in n, by dT]). Together with pi|) it implies that, for suffi- 
ciently large x and n > N, 

n-N n-N 

P{Wn > X} > {l-2e)Y E ^{^nkl}- 

k=l l=k+l 

Letting now n — t- oo, we derive from Lemma [5] the following lower bound, for all 
sufficiently large x: 

1 _ p poo 

P{W >x}> — — {Bi{x) f + b' / Bi{x + ya)B{x + y{a - b'))dy 

a{2a — b') L Jq 
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Note that, for any b' < b < a, 

foo a — b f°° 

/ Bi{x + ya)B[x + y{a-b'))dy> -/ Bi{x + ya)B[x + y{a -b))dy. 

Jo a - b' Jo 

We complete the proof of the Theorem by letting e J, 0. 

4. The maximal stability case: an upper bound 

Theorem 4. Assume b G (0,a). Suppose that the distribution Bj is subexpo- 
nential. Then, as x ^ oo, 



P{W >x}< '^'^^^ 



a{2a-b) 



{Bi{x)f + b / Bi{x+ya)B{x+y{a-b))dy 
Jo 



By Lemma [H it is sufficient to prove this upper bound for the queueing 
system D/GI/2 with deterministic input stream. So, let the interarrival times 
Tn be deterministic, i.e., r„ = a. Let an^ and an \ n > 1, be independent 
random variables with common distribution B. In this Section, define the service 
times an recursively. For that, we have to associate workloads with servers. Put 
Ui = C/1,2) = (0,0) and introduce the recursion 

= ([/„ + ea„cj„ - ia)+ (12) 

where = 1 if Un,i < Un,2 and a„ = 2 if C/„,i > C4,2- If Un,i = Un,2, then an 
takes values 1 and 2 with equal probabilities independently of everything else. 
Note that Wn = R{Un) a.s. for any n = 1, 2, .... 

Now we can define o"„ by induction. Lideed, ao is chosen at random from 
the set {1,2}. Put (Tq = o"o"°^ Then Ui is defined by recursion (fT2l) with 
n = 0. Assume that C/„ is defined for some n > 0. Then a„ is defined, too. 
Put (T„ = (t1""^ and determine Un+i by ([12]). 

Due to the symmetry, for any n, 

P{a„ = 1} = P{a„ = 2} = 1/2. (13) 

Consider two auxiliary D/GI /I queueing systems which work in parallel: at 
any time instant r„ = na, n = 1, 2, . . . , one customer arrives in the first queue 
and one in the second. Service times in queue i = 1, 2 are equal to a-h . Denote 
by Wn \ i = 1, 2, the waiting times in the ith queue and put Wn^ = Wn'' = 0. 
Since b < a, both queues are stable. Let W^^^ be a stationary waiting time in the 
zth queue. By monotonicity, with probability 1, 

W„ < min(VFii), WP) (14) 

for any n > 1. Hence, 

W < mm{W^^\w'^'^^}. (15) 
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Lemma 6. The waiting times {Wn^^} and {wji^^} are independent. 

Proof follows from the observation that the input (deterministic) stream and 
service times in the first queue do not depend on the input (also deterministic) 
stream and service times in the second one. 

Provided Bj is a subexponential distribution, 

P{W^W > x} ^Bi{x) as X ^ oo. (16) 

a — b 

Then Lemma [6] together with (jlSp implies the following simple upper bound: 

P{W>x} ^ 1 
lim sup — — --— < —77. 17 

Remarks. For a GI/GI/s queue with a < b and subexponential distribution 
Bj, similar arguments lead to 

P{W>x} ^ 1 
lim sup — < 



iBiix)y - {a -by 

Introduce the events, for k < n, 

A^^l = {4'^>x + in-kya-b)}, 
A^nl = Wk^>^ + in-kya-b)}. 

Lemma 7 (See also [3l Theorem 5]). Provided the distribution Bj is subexponen- 
tial, for any fixed N, 



limsupPliy^^) > ^' n ^nk] = o{Bi{x)) as x-^ 00. 

k=l 

Proof. For any 5 > 0, consider the disjoint events 

n-l 

Ci!^ = {{^k^ >x + {n-kya-b + 6)}nf] {af^ < x + (n - j){a - b)}} . 

Due to the Law of Large Numbers, there exists M > N such that 

P{t^«>x|cii)}>l-5 
for any n > M and k < n — M and, by the limit at ([7]), 

P{C2} > {l-6)P{a'i^ >x + {n-k){a-b + 6)}. 
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The events C^^"*, k < n — M, are disjoint, hence, 



n-M n-M 



k=l k=l 

n-1 

.(1) 



> (1 - 5f >x + k{a-b + 5))}. 

k=M 



The latter imphes the following lower bound: 

n-M 



liHQinf p{m/W > X, U Ci^2} > (1 - E ^(^ + A:(a - 6 + ,5)) 

(1-5)^ 



fc=l fc=Af 
\2 

--Biix) 



a — b + 6 

as X ^ oo. Since A^^^^ 5 C*^!;^ and since M > N and 5 > can be chosen 
arbitrarily, 

n-N 



lim^inf p{vfW > x, |J A^^l} > ^— ^S/(x) as x ^ oo. 



k=l 

Together with (jl6p . it implies the assertion of Lemma. 

Proof of Theorem U] continued. Estimate (jl4p and Lemma [6] imply 

71-Af n-Af n-N n-N 

i>{w^ > X, n 4^ u n 4?} < p{^i^^ > wi^^ > X, n 4^ u n 4?} 



A:=l Z=l A:=l i=l 

n-N. 

< ^ 

k=l 

n-N. 



/i — iv 

p{TyW>x, n 4^}P{H^P>x} 



+p{Ty« > x}p{h^(2) > p ^i2)y 



1=1 

Applying now Lemma [7] and relation p6|) . we conclude that, as x — t- oo, 

n-N n-N 

limsupP{M^„ > x, fl 42 U fl 4?} = o{(Bjix)f). 
""^^ ^ k=i 1=1 ' 

Since 



n-N n-N n-N n-N 

n 4'2u n ^1 = n K.'u4?) = u (^Sn4?). 

A:=l Z=l kl=l k,l=l 
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we obtain the equivalent relation, as x — >■ oo, 



n-N 



limsupP Wn>x, U (a^2 n 4?) = o{{Bi{x)f). (18) 



k.l=l 



Fix e > and put b' = b + e. For any n and k < I < n, define 

^2 = {4'^ >x+{l-k)a+{n- l){a - b')}, 
D^ni = {<^?^>^ + in-l){a-b')}, 
Dnkl = {<^k^ + >2x+{l- k)a + (n - Z)(2a - b')}. 
For any n and I < k < n, define 

Di^^ = {a'j}^>x + in-k){a-b')}, 

L>2^ = {cTp^ > a; + (fc - Oa + (n - fc)(a - 6')}, 

Dnkl = Wk^ + > 2x + (A: - l)a + {n - k){2a - b')}. 



Denote 



We can derive an upper bound on the probability of the event {Wn > x} as 
follows: 

P{W„ > x} 



n-N n-N n-N 

< P{W„ >X, \J Fnkl} + P{Wn >X, \J FnkU U I 
k,l=l k,l=l k,l=l 



n-N 



+p{w^>x, uK^4f)} 



k,l=l 

= Pnl + Pn2 + Pn3- (19) 

Here the first term is not greater than 

n— 1 n— 1 n— 1 

Pnl < P{Wn >X, \J Fnkl} + P{Wn > X, [J F^kl] + P{Wn > X, \J F^kk} 

k.l = l k,l=l k=l 

k<l k>l 
= Pall + Pnl2 + Pnl3- (20) 

The third probability is negligible in the sense that 

n— 1 n— 1 

Pnis < p{ U (^2 n DSi)} < E p{^2}p{^S} 

k=l k=l 
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oo 



< B{x) B{x + k{a-h- e)) 

k=l 

< B{x)Bi{x)/{a -b-e) = o{{Bi{x)f) (21) 



as X — )• oo, since B{x) = o{Bj{x)). The first probability in ([2U|) admits the 
following upper bound: 

n— 1 n— 1 

fc=l «=fc+l 

n— 1 n— 1 

k=l l=k+l 

For Si, we have the following inequality and equalities: 

n-l 



k<l 

n—1 2 'I— 1 

= J2 = l}P{^S.^i?'^n«} = 2 E nPnkl}, (22) 

fc,i = l fe,! = l 

fc<i k<l 

by independence of the event {a^ = 1} from D^^j^, D^^ and -D^fci and by the 
symmetry (jl3p . The sum S2 is not greater than 



fc=i 

71—1 n — 1 

E P{D^^hP{Wn >x,ak = 2}< P{Wn > x} E P{^S}- 



k=l k=l 

Hence, S2 = o(P{Wn > x}) as x — t- 00 uniformly in n. Combining the latter fact 
with estimate ([22]) for Si, we get 

-j^ n—1 

i^nll < 2 E ^{^nkl} + oiP{Wn > x}). (23) 

fc,i=l 
k<l 

Taking into account the equality P„ii = -Prai2, we obtain from (j20]) . ([2T]) and ([23]) 
the following estimate: 

n-l 



Pnl < E P{i^nH}+o((57(x))^ 



fc<« 
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as X — )• oo uniformly in n. Now applying the calculations of Section 3.3 we can 
write down the following estimate, as x — >• oo: 

IimsupP„i < \ . {Bj{x) f + b' / Bj{x + ya')B{x + y{a - b'))dy . 

n^oo a{2a — b') L Jo J 

(24) 

It is proved in (|18p that, uniformly in n, 

P„3 = o{(Bi{x)f) as X ^ CO. (25) 

We have 



n-N n-N n-N 



k,l=l k,l=l k,l=l 

Thus, 

n~N 

Pn2 < E nwn > X 1 4\\4?,F„,,}p{4'2 n4¥}. (26) 

k,l=l 

Conditioning on Wnk and Wni yields, for any tt; > 0, 

P{Wn > X I 4\\A2\F„fci} < P{Wn >x\Wkl< W,Wi2 < U;, ^fc , , ^nfc J 

+P{Wkl >W}+ P{Wl2 > w}. 

Since b < 2a, the two-server queue is stable and, in particular, the sequence of 
distributions of random variables (Wni,Wn2) is tight. It means that, for any 
fixed e > 0, there exists w such that, for any A; > and / > 0, 

P{Wfci >w}<e and P{Wi2 > w} < e. 

Also, from the stability and from Corollary [H for any fixed e > and w > 0, 
there exists such that, for any n > N and k, I < n — N, 

P{Wn >x\Wkl< W,Wl2 < < ^- 



Combining these estimates we obtain from ([26]) . 

n—N n—N 

Pn2 < 3. E nA^nl n 4?} = 3< E P{4^ })'■ 

k,l=l k=l 

Hence, 

Pn2 < 3e{j2B{x + kia-b'))) < ^-^(Bjix)f. (27) 

Since the choice of e > is arbitrary, relations (fT9l) ~ ([25]) and (|27|) imply the 
conclusion of Theorem HI 
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5. The minimal stability case: lower bounds 

Theorem 5. Let h G (a, 2a) and the integrated tail distribution Bj be long 
tailed. Then the tail of the stationary waiting time satisfies the following inequal- 
ity, for any fixed 6 > 0: 

^rr.r . 1 + o(l)— i'b + 6 \ 

P{W >x} > -t-BA- x] as x-^oo. 

^ ^ - 2a-b \b-a ) 

Notice that if 6 G (a, 2a) then -^>2. 

Remark 4. By use of similar arguments, one can get the following result 
for an s-server queue, s > 2: if the integrated distribution Bj is long tailed and 
6 G ((s — l)a, sa), then, for any 5 > 0, 

sa — b V — (s — \)a 
Theorem [5] implies the following 
Corollary 6. Assume that Bi G ,y.'%'f. Then, as x ^ oo. 

In the case b G [a, 2a), one can also derive a lower bound which is similar to 
([6]). More precisely, assume b G [a, 2a). Then introduce another two-server queue 
with the same service times and with inter-arrival times = ct„, where c > b/a. 
For this queue, denote by W a stationary waiting time of a typical customer. 
Due to monotonicity, P{1^' > x} < P{W > x} for all x. Applying Theorem [3] 
and Remark [H we get the following lower bound for the case b G [a, 2a): if the 
integrated tail distribution Bj is long-tailed, then, for any c> b/a, 

P{H'>.)>(l + „(l))^j5|2i±^(B,W)^. (28) 

Proof of Theorem \5\ By Lemma it is sufficient to prove the lower bound 
for the queueing system D/GI/2 with deterministic input stream. Let the inter- 
arrival times Tn be deterministic, i.e., r„ = a. For any 5 > 0, set e = . Put 
b' = b — e and = For any /c G [l,n — N], consider the events 

Ank = Wk >2x + (2a - b'){n - k)}, 

n 

Cnk = f]Wi<2x + (2a - b'){n - I)}. 
1=1 



as x 



oo. 
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Since Ecj is finite, 

oo 

F{Cnk} < Pi^i > + (2a - b')l} = 0(Bi{2x)) (29) 
1=1 

as X — )• 00 uniformly in n > 1 and k < n. Since tlie events A^k H C^^, A; G [1, n], 
are disjoint, we obtain 

n-N 

P{Wn >X}> P{Wn > X,Ank,Cnk} 

k=l 

n-N n-N 

> J2 > ^'^nj - E nAnk,Cnk}. 

k=l k=l 

Since the events Ank and Cnk are independent, we get 

n-N n-N 
F{Wn >X}>Y > X,Ank} " SUpP{C„fc} J] P{A„fc} 

fc=l fc=l 

= J2 > ^ I Ar,k}P{Ank} - 0(1) 5] P{A,fc} (30) 

fe=l fc=l 

as X —7- 00 uniformly in n > 1, by (j29p . The event implies the event 

> 23; + (2a - 6')(n - A:) - a. 
Thus, it follows from Corollary [5] that 

P{Wn > X I Ank} ^ 1 

as x — )• 00 uniformly in n and k < n — N . Therefore, we can derive from (j30p the 
estimate 

n-N 

P{W >x}= lim P{M^„ >x} > (1-e) lim V P{A„fe} 

k=l 

00 

= (1 - e) E S(2x + (2a - b')k), 

k=N 

which is valid for all sufficiently large x. Since the tail Bi{v) is long-tailed, 

00 -1 

y 'B{2x + (2a - h')k) 77^/ (2x + (2a - b')N) 

tN 2a -6' 



= 2^^K6^^J = 2^^K6^"J 
as a; — )• 00. The proof is complete. 
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6. The minimal stability case: an upper bound 

Theorem 6. Assume b G [a, 2a). Let both B and Bj be subexponential dis- 
tributions. Then the tail of the stationary waiting time satisfies the following 
inequality, as x ^ oo: 

p{w >x} < V-37r^^(2x). 

Remark 5. By use of the same arguments, one can get the following result 
for any s-server queue, s > 2: ii Bj £ y and b < sa, then 

P{W > x} < ^ ^^]^ Bi(sx) as X oo 
sa — b 

provided that either (i) ai > (s — l)a a.s., or (ii) B G 5^ . 

Remark 6. For an s-server queue, Foss and Chernova [10] have proposed 
another way of obtaining upper bounds; it is based on comparison with a queue 
with the so-called cyclic service discipline. 

Proof of Theorem [6l From Lemma O it is sufficient to consider the case of 
constant interarrival times r„ = a only. Put M„^o = and 

Since b > a, Mo,„ — )■ oo a.s. as n — )• oo and, due to the Law of Large Numbers, 

Mo,„ 



n 



a a.s. (31) 



and in mean. Note that EMg^n > n{b — a), since Mo^„ > ctq + . . . + — na. 
For any given e > 0, choose an integer L > such that 

€[b-a,b-a + e). (32) 

Consider any initial workload vector Wq = (Wo,i,Wo,2) > 0. Put Z„ = 
Wn,i + Wn,2- Since the increments of the minimal coordinate of the waiting time 
vector is not greater than the increments of Mo^n, 

Wi^n — Wifl < Mo,n for any n. 

Hence, provided Wn^2 ^ we have the inequality 

Zn+l — Zn < Mo,n+l — Afo,n — «• 

If Zq > 2aL, then Wo,2 ^ clL and, for n = 0, . . . , L — 1, Wn,2 ^ o,{L — n) > a. 
Therefore, if Zq > 2aL, then 

Zl< Zo + Mo,L - aL. 
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Monotonicity implies, for any initial vector Wo, 

Zl < max{2aL, Zq} + Mq^l — o-L. 

Thus, the following inequalities are valid for any n: 

Z{n+i)L < max{2aL, Zni} + MnL,L - aL. (33) 

Consider a single-server queue with i.i.d. service times an = ^nL,L and con- 
stant inter-arrival times r„ = La and denote by Wn a waiting time of nth cus- 
tomer. This queue is stable since b = E^i < aL = a. Put Wq = 0. Assuming 
that Zq = 0, we can derive from (j33p the following bounds: for all n = 0, 1, . . . , 



ZnL < 2aL + Wn a.s. (34) 

Denote G{x) = P{ao > x}. We show that integrated tail distribution Gj is 
subexponential one. We need to consider only the case L > 1. Note first that 

(To + . . . + aL-i - La < ao < (To + . . . + ctl-i a.s. (35) 

Since the distribution of ai is assumed to be subexponential, the asymptotics for 
the lower and upper bounds in the latter inequalities are the same: as x — )• oo, 

L-l L-l 

fji - La > x} ~ CTi > x} ~ LS(x). (36) 



Therefore, the tail G{x) has the same asymptotics and G is a subexponential 
distribution. Thus, 

Gi{x) = / G{y)dyr^ LBiix). (37) 

and Gj is a subexponential distribution, too. Thus, by classic result ([T|) for the 
single server queue, the steady state distribution of the waiting time Wn satisfies 
the following relations, as x — )• oo: 

lim P{Wn > x} K.Gj(x) < 1 — rrGiix) ~ —4 — Bi(x), (38) 

n^oo ^ ^ a-b {2a-b-e)L ^ ' 2a-b-e ^' ^ ^ 



by ([32]) and (laiD. Since Z„ = Wn,i + Wn,2 > 2W„,i, 

P{W > x} = P{2W > 2x} < lim P{Z„l > 2x}. 

n— i>oo 

Now it follows from ([SD and ([38]) that 

P{W > x| < lim P{Wn >2x- 2aL} 

since Bj is long-tailed. Letting e J, concludes the proof. 
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7. The minimal stability case: exact asymptotics 

In this Section, we prove Theorem [2l First note that, as fohows from (I28j) . the 
tail P{1^ > x} may be heavier than that in Theorem [21 in general. For instance, 
this happens if 

B[(- x) = o(b'^t(x)) as X — )• oo. (39) 

vo — a ^ 

Assume b E (a, 2a) and consider, for example, a service time distribution with the 
Weibull integrated tail Bi{x) = e-^^ /3 e (0, 1). Then 1^ holds if {-^f > 2. 

Proof of Theorem [2l Since Bj £ ^^'f, both the lower bound in Theorem 
[5] and the upper bound in Theorem [6] are of the same order, 

Bii2x) = 0(Bi[-^x)). (40) 

We use the notation from the previous Section. In particular, we fix e > 
and choose L satisfying (j32p . For any constant c > 0, ()35p implies 

L-i , L-i . 

U WkL+i > x + La + {L- i)c} C J ^ o-fc^+j - La > x I C {5=^ > x}. 

i=0 ^ i=0 ^ 

Therefore, from ([35]) and ([36]) . 

P{{5fc > a:} \ U {akL+i >x + La + {L- i)c}} = o(B{x)). (41) 

2=0 

Take c = (a-b)/L. By ([Ml), 

P{VF >x} = lim P{1^„L 1 > 2;} = lim P{W„l 1 > x, > 2x - 2aL}. 

n— >oo ' n— s>oo ' 

Standard arguments concerning how large deviations in the single server queue 
Wn occur imply the relation 

n-l 



P{W > x} = lim J2 P{^nL,i > x,ak>2x + {n- k){a - h)} + o(S/(2x)) 

fc=0 
nL-l 

lim y P{VF„L,i > x,Oi>2x + {n- i)c] + o(^/(2x)). 



n—KX) 

1=0 



by dH]). Now it follows from ([32]) that 

n.L-l 

P{W > x} < lim y PlWnL 1 > x,o-i > 2x+(n-ii)(2a-6-e)} + o(B7(2x)) 

2=0 
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nL 

< ^}}^Y.^i^nL,i > x,anL-j > 2x+j{2a-b+e)} + sO(Bi{2x)) + o(Bj{2x)) 

/Nil-e) nL \ 

= lim V + V ]+eO{Bi{2x))= \\m{ll^ + i:2)+eO(Bi{2x)), 

\ j=l j=Ar(l-e)/ 

where N = x/{h — a). The second term admits the followmg esthnate 

oo 

S2 < ^{^ >2x + j{2a-b)} 

j=N(l-e) 

~ — ^—Bi{2x + N{l-£){2a-h)) = — ^—bA-^x - e ^" ~ ^ x\ . 
2a-b ^ ^ " 2a-b - a b-a J 

It follows from Bj G J^^'f that, for any S > 0, there exists £ > such that 

which coincides with the lower bound in Theorem [5l 

Now consider the first term Si. Since the queue is stable, one can choose 
K >0 such that P{W„,2 < -fC} > 1 - e for ah k. Then 

Nils) 

Si < ^ P{WnL-j,2 > K, anL-j > 2x + (2a - 6 + e)j} 
i=i 

Nils) 

+ J2 ^{WnL,l > X, WnL-j,2 < K, UnL-j > 2x + (2a - 6 + e)j} 

i=i 

= Si^i + Ili^2- 

We have 

Nils) 

Si,i = J2 ^{WnL-j,2 > K}P{ai > 2x + (2a - 6 + e)j} 
i=i 

oo 

<eY, ^Wi > 2x + (2a - b)j} < ^^Bj{2x). 
Note that if WnL-j,2 < K, then WnL,i < K + MnL-j+i,j-i- Therefore, 

Nils) 

Sl,2 < J2 ^WnL-j >2X+ (2a - b)j, K + MnL-j+l,j-l > x} 

i=i 

Nils) 

= ^WnL-j >2x + (2a - b)j}P{K + Moj_i > x}. 



26 



S. FOSS AND D. KORSHUNOV 



Since the sequence Mq j stochastically increases, 



Si,2 < P{K + Mo,^(i_,) >x}J2 P{^i > 2x + (2a - b)j} 

1 - 

2a - b' 



< P{Mo,^(i_,) >x- K},^—^Bi{2x). 



Since 



A^(l - e) 1 - e 
we have by ([3T]) 



X — K b — a 

> b — a as X — )• oo, 



Thus, we have shown that the upper bound for P{VF > x} is not bigger than 
the lower bound in Theorem [5] plus a term of order 

/ / bx 
{e + 6)0{Bi{2x)) < {e + S)0(B/ 



b — a 



due to (PUI) . Since e > and 6 > may be chosen as small we please, the proof 
of Theorem [2] is complete. 

8. Tail asymptotics for the two-dimensional w^orkload vector 

Denote by = {W^, W2) a weak limit for the vectors Wn as n — t- c«. Clearly, 
W = W^. 

[8].l. Maximal stability case. First, we obtain simple lower and upper 
bounds which are equivalent up to some constant. Second, we give (without a 
proof) a result related to the exact asymptotics. 

Theorem 7. Let b < a and Bj £ ^ . Then, as x, y ^ 00, x < y, 

P{< > x,W^ >y}> ^d:^Bj{x)Bj{y). 
//, in addition, Bj £ ^ , then 

P{< > x,W^ >y}< l±^Bj{x)Bj{y). 

Proof. Fix e > and put a' = a + e. For k, I < n, k ^ I, define the events 
Anki and Cnki by the equalities 

Anki = {o-fc > X + {n- k)a, ai > y + {n - Z)a'| 
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and 



7t 

Cnki = {~\ [cFj <x + {n- j)a'y 



Note that the events Arifci ^ C„fc; are disjoint for different pairs (fe, I) and 

n n 

P{Ty„l > X,Wn2 >y} >Y^ Yl P{^nl > X^Wn2 > V , A^kl , Cnkl} ■ 

k=ll=k+l 

Then the same calculations as in Subsection [3l3 imply the estimate, as x, y — )• oo, 

n—l n— 1 

P{Wni > X, Wn2 > y} > (1 + 0(1)) ^ J] B{x + {n - k)a')B{y + (n - l)a') 

k=i i=k+i 

n— 1 n—k—1 

= {l + o{l))J2 B{x + ka')B{y + la'). 

k=l 1=1 

Hence, 

oo oo 

P{< > X, > y} > (1 + 0(1)) E E^(^ + ^"'My + ^«') ~ ^/(x)5Ky)/. 

k=l 1=1 



12 

a 



Proceed to the upper bound. Due to construction of the majorant {Wn'^ , Wn'' ] 



and the lower bound is proved. 

Proceed to the upper bound 
in Section m we have the inequality 

P{< > X, >y]< lim [P{W^^i) > X, VFi^) > y] + P{T^^i) > y, vr^2) > ^1 

n— >-oo L 

= 2 lim P{W(^) > x}P{W(2) > y}. 

n— ^oo 

Together with (jl6p it implies the desired upper bound. Theorem [7] is proved. 

Turn now to the exact asymptotics. Below is the result. The proof is rather 
complicated and will be presented in another paper. Denote 

/-oo 

R{x, y) = Bi{x)Bi{y) + h Bj{y + za)B{x + x{a - h))dz. 

Jo 

Recall that Theorem 1 states that P{W^ > x} ~ R{x, x)/a{2a-h) given Bi G y. 
Theorem 8. Assume b < a and Bj € o5^. Let x,y ^ oo, x <y. Then 

F{W^>x,W^>y} ~ ^ R{y,y) + ^{R{x,y)-R{y,y)). 

a{2a — 0) 
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[8].2. Minimal stability case. We prove the following 

Theorem 9. Assume a < b < 2a, B £ Y, and Bj G J^^f. Let x,y — )■ oo 
in such a way that y /x ^ c G [l,oo]. Then 

Proof. Start with the case c = oo. From Theorem 10 in [3], one can get the 
following: 

Corollary 7. Assume b G (a, 2a). If B G y and Bi £ .5^ , then, as y ^ oo, 

P{< >y}- -Bi{y) + -^^Bj(y-±-). 

a a[2a — b) ^ b — a' 

It is clear that 

V{Wf > X, >y}< P{W^ > y}. 
On the other hand, for any = 1, 2, . . . , 

P{< > X, >y} = lim P{Wn,i > X, Wn,2 > y} 

n— ^oo 

n-1 

> lim F{Wn-N.2 >y + Na, V (dj - tj) > x] 

j=n-N 
N 

= lwa^P{Wn-N,2 >y + iVa}P{^(c7j - Tj) > x}. 

j=l 

Fix e > 0. Put iV = N{x) = x{l + e)/{b - a). Then by LLN 

N 

for all sufficiently large x and, as n — )• oo, 

P{Wn-N,2 >y + Na}^ P{W^ >y + Na}. 

Since Bi G ^^Y, 

F{W^ >y + Na} ~ P{W^ > y} as y ^ oo. 

By letting e — )• 0, we get the result. 

Now consider the case c < oo. If c = 1, then the result follows from Theorem 
[2j Let c G (l,oo). We give here only a sketch of the proof, by making links to 
the proof of Theorem [2j 
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Since 

P{W^ >y}< P{W^ > X, >y}< P{W^ > x} 

and 

P{W^ > y} ~ P{W^ > cx} >{K + o{l))P{W^ > x} 

where K = mftBj(ct)/Bj{t) > 0, one can get from the proof of Theorem [2] 
the following equivalences: for Nx = x/{b — a), Ny = y/(b — a), and for e G 
(0,1-1/VH), 

P{W? > x,W^ > y} 

n-N^{l-e) 

lim > X, Wn,2 >y,(Ji>2x + {n-i){2a-b)} + eO(5/(2x)) 



n—^oo 

i=l 



/n-Ny{l+e) n-Af:c(l-e) \ 

= lim V + V ] +eO{Bi{2x)) = {1:1 + ^2) +eO{Bi{2x)). 

\ i=l i=n-Ny{l+e)/ 

Choose K >Q such that P{Wn,2 > K} < e for ah n. Then 

n-iVa;(l-e) 

S2 = lim V P{W^,2<K,Wn,l>x,Wn,2>y,(yi>2x + {n-i){2a-b)} 

i=n-Ny(l+e) 

+eO(Bi{2x)). 

From Lemma 2 and its Corollaries, 

Ny{l+e) 

S2 = (l + o(l)) PWi>y + ja} + eO(Bj{x)) 

j=N4l-e) 



due to -B/ G J^^'V. From Lemma 3 and its Corollaries, one can also conclude 
that, for ? < n - Ny{l + e), if < 2y + (n - i){2a - b - e) and Wi^2 < -f^, 
then, with probability close to one, both coordinates of the vector W„,2) 
take values less then y for all sufficiently large n. From the other side, if cij > 
2y + (n — i){2a — b + s), then, with probability close to one, y < Wn,i < Wn,2- 
Therefore, 

00 

El = (1 + 0(1)) P{^i >2y + i(2a-6)} + eO(;BKx)) 

j=Ny{l+e) 

-'^"^'hA^)+eOmx)). 



2a — b \b — a' 
Summing up the terms and letting e and 5 — )• concludes the proof. 
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9. Comments on stationary queue length 

Let Qn be a queue length viewed by an arriving customer n, and Q its stationary 
version in discrete time (i.e. Palm-stationary). Due to the distributional Little's 
law, 

P{Q > n} = Y{W > Tn] 

where W is the stationary waiting time, T„ = ri + . . . + r„, and W and do not 
depend on each other. When a distribution of W is long-tailed, the asymptotics 
for P{VF > Tri}, n — )• oo, have been found in [2] and in [TT]. If, in addition, t„ 
has a non-lattice distribution, there exists a stationary distribution G for a queue 
length in continuous time. Then, from Lemma 1 in [TT], 

G{n) ~ P{<5 > n} as n — )■ oo. 
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